Higher order proton structure corrections to the Lamb shift in muonic hydrogen 



Carl E. Carlson 1 ' 2 and Marc Vanderhaeghen 3 

1 Helmholtz Institut Mainz, Johannes Gutenberg-Universitat, D-55099 Mainz, Germany 
z Department of Physics, College of William and Mary, Williamsburg, VA 23187, USA 
3 Institut fiir Kernphysik, Johannes Gutenberg-Universitat, D-55099 Mainz, Germany 

(Dated: January 31, 2011) 

The recent conundrum with the proton charge radius inspires reconsideration of the corrections that 
enter into determinations of the proton size. We study the two-photon proton-structure corrections, 
with special consideration of the non-pole subtraction term in the dispersion relation, and using fits 
to modern data to evaluate the energy contributions. We find that individual contributions change 
more than the total, and present results with error estimates. 
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The recent measurement of the proton charge radius 
using the Lamb shift in muonic hydrogen [1] has given a 
value that is a startling 4%, or 5 of the previous standard 
deviations, lower than the values obtained from energy 
level shifts in electronic hydrogen 13] or from electron- 
proton scattering experiments (j- \m- Specifically, the 
new muonic hydrogen measurement Jlj] gives 



The (9 (a 5 ) term was given by Friar ^ as, 



R E = 0.84184 (67) fm, 
compared to the COD ATA value 0] 

R E = 0.8768 (69) fm, 
or the latest electron scattering value 0] 
Re = 0.879 (8) fm, 
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where we have added in quadrature the several uncer- 
tainties given in Q]. 

The promise of the muonic hydrogen measurement 
was that, because a muon would orbit closer to the 
proton than an electron, the effect of the proton struc- 
ture on the energy level splittings would be enhanced 
and a more accurate proton radius could be obtained. 
Based on the quoted error limit, that promise has been 
achieved. However, the discrepancy from the previous 
results requires an explanation, and invites a reconsid- 
eration of the theoretical corrections that are involved in 
connecting the experimental energy shift to the proton 
charge radius l3,|6[]. In this note, we will focus on one of 
the corrections, namely the order a 5 proton size correc- 
tions to the Lamb shift. 

The leading O (a 4 ) and O (a 5 ) proton structure contri- 
butions to the hydrogen Lamb shift are often given as 
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where </> 2 (0) is the square of the nS-state wave function 
at the origin (which contains a factor a 3 ) and m r is the 
lepton-proton reduced mass. The quadratic term was 
obtained non-relativistically in |@], and one can verify 
from a relativistic calculation that Re is indeed the pro- 
ton charge radius JjJ. 
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where Pe is the charge density of the proton itself. Friar 
called R 3 2 j the third Zemach moment, because it is rem- 



iniscent of an integral found by Zemach MlOU in the re- 
lated context of hyperfine splitting. 

In modern times, one should calculate the 0(a 5 ) cor- 
rections field theoretically using the diagram shown in 
Fig. as has been done by Pachucki JTlLFl^h and by oth- 
ers 111 311 . We wish to reexamine the calculation here, for 
the purpose of better assessing the connection between 
the elastic and inelastic contributions, and to better eval- 
uate the subtraction term needed in a dispersion relation 
that is part of the work. 

The calculation of the elastic and inelastic contribu- 
tions should be done together. Perhaps in the future 
a direct QCD calculation will be possible, and there is 
an exploration of the hadronic corrections to the Lamb 
shift using crural perturbation theory fl3l , but for the 
present to obtain the required accuracy the calculation 
needs to be done dispersively, connecting the off-shell 
Compton scattering which is the hadronic side of the 
diagram to information obtained from electron-proton 
scattering. In particular, done that way the elastic con- 
tributions require no (non-existent) knowledge of form 
factors for situations where a proton is off-shell. It also 
means that certain non-pole contributions to the Comp- 
ton amplitudes are not picked up by the dispersive cal- 
culation and in fact do not contribute. 

We also analyze more concretely the subtraction func- 
tion that appears because one of the dispersion rela- 
tions does not converge if unsubtracted. The subtrac- 
tion function depends on the photon four-momentum 
squared, Q 2 , and its value at Q 2 = is given in terms of 
the proton magnetic polarizability. Its Q 2 dependence 
can be estimated by calculating a two-pion loop contri- 
bution which couples to the nucleon as a scalar. One 
does not need to use a Q 2 dependence assumed given 
by the nucleon electromagnetic form factor, as has been 
done previously. 
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FIG. 1: The box diagram for the 0(k 5 ) corrections. 



The Feynman diagram for the two-photon proton- 
structure correction to the Lamb shift is shown in Fig. [T] 
To the level of accuracy needed here, all external lines 
have zero three-momentum. The blob corresponds to 
off-shell forward Compton scattering, given in terms of 
the Compton tensor 
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where q 2 = —Q 2 , v = p ■ q/M, and M is the nucleon 
mass. A spin average is implied and the state normal- 
ization is (p\p') = (27r) 3 2E<5 3 (p — p'). The functions 
Ti, 2 {v,q 2 ) are each even in v and their imaginary parts 
are related to the structure functions measured in elec- 
tron or muon scattering by 

ImT 1 (v,Q 2 ) = ^F 1 (v,Q 2 ), 

ImT 2 (v,Q 2 ) = i ; F 2 (t/,Q 2 ), (7) 

with v > and where F\ 2 are standard 

After doing a Wick rotation, where ^0 = 'Qo an d Q = 
q, one obtains the 0(a 5 ) energy shift as 

AE = ^(0)/A2 

(Q 2 + 2Q 2 )T 1 (zQ , Q 2 ) - (Q 2 - Q 2 )T 2 (iQ , Q 2 ) 



Q4 (Q 4 +4m 2 Q 2 ) 
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where m is the lepton mass, and <p 2 {Q) = m 3 a 3 / (7m 3 ) 
with m r = mM/(M + m). 

The T; are obtained using dispersion relations. Regge 
arguments 1 16| suggest that T 2 satisfies an unsubtracted 
dispersion relation in v at fixed Q 2 , but that T\ will re- 
quire one subtraction. Before proceeding, we will note 
that the Born terms, obtained from the elastic box and 
crossed box of Fig. [2] and the vertex function TP = 



l fl h(Q 2 ) + (i/2M)a' iV q v F 2 (Q 2 ) for an incoming pho- 
ton, are 

T?(qo,Q 2 ) 
T!(qo,Q 2 ) 



Q'GUQ 2 ) 



AtzM \ (Q 2 - ie) 2 - AM 2 q 2 

MQ 2 G 2 (Q 2 ) + t p G 2 m (Q 2 ) 
n(l + T p ) (Q2 - ie) 2 - 4M 2 q 2 



F 2 (Q 2 ) 
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where T« = Q 2 / (4M 2 ), and the electric and magnetic 
form factors are 

G E (Q 2 ) = F 1 (Q 2 ^ 



G M (Q 2 



h(Q 2 )+F 2 (Q 2 ). 
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The Born terms are reliable for obtaining the imaginary 
parts of the nucleon pole terms, but not reliable in gen- 
eral, since the given vertex assumes the incoming and 
outgoing nucleons are both on shell. 

Calling the first term in Tf the pole term, one can split 
the whole of T\ into pole term and non-pole terms, 
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The pole term alone evidently allows an unsubtracted 
dispersion relation, and this term calculated from the 
dispersion relation simply reproduces itself. With a once 
subtracted dispersion relation for T\, one has 



Ti(q ,Q 2 ) = T : 



f o VQ 2 )+Ti(0,Q 2 ) 

r dv Uv,Q 2 ) 
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The nucleon pole is isolated in T? e a nd the integral 
begins at the inelastic threshold i/y, = (2Mm n + m 2 + 
Q 2 )/(2M). Similarly, as contains only a pole term, 



T 2 (qo,Q 2 ) = T 2 B (q ,Q 2 ) + =- 



2/T 



dv 



F 2 (v,Q 2 ) 
v 2 -q\ 



With 



we obtain 



A£ 



subt 



AE 



Ana. 



AE suht + AE ind 



AE e 



(13) 



(14) 



, <K(0), 
tn Jo 



dQ 2 7ifo)~ 



Q 2 



Tt(0,Q z ), (15) 



k k 



P V 




FIG. 2: Elastic contributions to the box diagram. 
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where t = v 2 /Q 2 and = Q 2 /(4wz 2 ). The auxiliary 
functions are 

7i(t) = (1-2t)((1 + t) 1 / 2 -t 1 / 2 )+t 1 / 2 , 

72 (t) = (1 + t) 3 / 2 -t 3 / 2 -^/ 2 . (18) 

Both are monotonically falling functions, reducing to 1 
at t = and falling like t _1//2 at large t. Also 

1 
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The subtraction function T(0, Q 2 ) has unphysical ar- 
guments, excepting the point Q 2 = 0. It comes from the 
excitation of the proton, and can at low Q 2 (and low v, 
in general) be described using the electric (a^) and mag- 
netic (/3m) polarizabilities and the effective Hamiltonian 



n = -Una E E 2 -hnl5 M B : 



(20) 



For small v and Q, this gives 

— v 2 O 2 

lim Ti(v,Q 2 ) = — ( aE + /3 M ) + i^/3 M . (21) 

v 2 ,Q 2 ^o e e 

The v 2 term is shown to connect to known results in 
another context fl^|, and the Q 2 term was obtained by 
Pachucki Jll]] . With the above result, the integral over 
T\ (0, Q 2 ) converges at the lower limit. 

For higher Q 2 , the subtraction function comes from 
non-nucleon-pole contributions, and the forward am- 
plitude is dominated by low mass intermediate states. 
With the Q 2 — » limit fixed in terms of /3m, we estimate 
the Q 2 dependence from pion loop contributions where 
the two-pion state has a scalar coupling to the nucleon, 
as illustrated in Fig. [3] 

With standard Feynman rules for scalar QED and an 
effective gsNcp^tpN coupling for the lower vertex, one 
obtains from these terms 
gs 
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FIG. 3: Diagrams used for estimating the Q 2 dependence of 
the non-pole part of the subtraction term. 



where, with A = 4m 2 / Q 2 , 
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We can identify /3 m = ags/487r 2 m 2 , and obtain 
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Ti(0,Q 2 ) = -^Q 2 F loop (Q 2 ). 
The Particle Data Group gives jl5ll , 

^M = (1-9 ±0.5) x 10~ 4 fm 3 . 
However, according to some recent analyses, 



/3m 



(4.0 ± 0.7) x 10~ 4 fm 3 [18] 
(3.4 ±1.2) x 10~ 4 fm 3 [19,20] 



(24) 
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Using the subtraction function from Eq. I I24I I, we find 



AE subt = 5.3 }teV x 



(3.4 x 10- 4 fm 3 ) 



(27) 



Much of the support for the integral is at low Q 2 , being 
controlled by 7 i as well as by the Q 2 dependence from 
the pion loop, and half the contributions to AE subt come 
from Q 2 < 0.04 GeV 2 , albeit with a long tail. 

Refs. Jll and (H found AE s " bt to be 1.8 and 2.3 pieV, 
using (6m = 1.5 and 1.9 x 10~ 4 fm 3 , respectively, and 
using a Q 2 falloff related to the nucleon electromagnetic 
form factor. For the same /3m, our results are about 30% 
larger due to having flatter Q 2 falloff. 

One can also consider inserting a form factor F n for 

each incoming photon coupling to pions, modifying the 

subtraction function of Eq. l[24j by multiplying it with 

F n (Q 2 ) 2 . Obtaining F n from the fit of |22], we find 
AE subt = 3 8 ?(eV 

The inelastic contributions depend on F(i,2) { v i Q 2 )/ 
and good data in the low-Q 2 and resonance region is 
available from Jefferson Lab. Analytic representations 
of this data are given by Christy and Bosted |23||. in a 
fit valid for < Q 2 < 8 GeV 2 and W from threshold to 
3.1 GeV, where W is the final hadronic mass for inelastic 
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ep scattering, W 2 = M 2 + 2Mv - Q 2 . From the Bosted- 
Christy region, we obtain a —12.2 ^eV contribution to 
AE" ul . We also use the fit of Capella et al. HI, valid for 
data at low and intermediate Q 2 above the resonance 
region, specifically < Q 2 < 5 GeV 2 and W > 2.5 
GeV. This gives a —0.5 }ieV contribution using l24Tl for 
W > 3.1 GeV in the allowed Q 2 region. Contributions 
from higher Q 2 are quite small (on the order of 0.002 }ieV 
from Q > 5 GeV 2 and W > 3.1 GeV). We thus have 

k E inel = -12.7 jieW . (28) 

Refs. 03] and fH quoted -13.9 and -13.8 jieV for this 
contribution. 

The elastic contribution depends on the nucleon form 
factors, for a selection of modern form factors we get 

C -27.8 jieV Kelly [25] 
AE el = \ -29.5 ^eV AMT [26] . (29) 

[ -30.8 ^eV Mainz 2010 [4,27] 

Ref. jli n q uoted —23 pieV using the Simon et al. form 
factors | 28|| from 1980. However, the main difference be- 
tween our results is not due to the newness of the form 
factors, but rather to our exclusion of the non-pole con- 
tributions from the elastic contributions. The non-pole 
contributions would be a positive 4.7 (Kelly) or 4.8 jieV 
(AMT or Mainz 2010) contribution were they included. 

Table U summarizes our numerical results (selecting 
the AMT form factors |26|| for the elastic terms) and com- 
pares them to earlier work. (Ref. l2lll did not calculate 
the elastic term, so we carried over the result from jllll .) 

TABLE I: Numerical results for the O (a 5 ) proton structure cor- 
rections to the Lamb shift in muonic hydrogen. Energies are in 



peV. 


(jieV) 


this work 


Ref. [1L 12] 


Ref. [21] 


AE subt 


5.3 ±1.9 


1.8 


2.3 


AE' ne < 


-12.7 ±0.5 


-13.9 


-13.8 


AE el 


-29.5 ±1.3 


-23.0 


-23.0 


AE 


-36.9 ± 2.4 


-35.1 


-34.5 



Regarding the uncertainties, for the subtraction term 
energy, we took fi^ from l2(ill , and propagated their er- 
ror limits, which are large enough to accommodate the 
other two /3 m values. The inelastic energy, comes mainly 
from |23|l . which states that most of the data points are fit 
to within 3%. The data itself typically had 3% error lim- 
its, and we added these two errors in quadrature. For 
the elastic term, we estimated the error from the spread 
between the two newer form factor fits that we used. We 
added the errors in quadrature to obtain the total error. 

Our results are similar to previous results in aggre- 
gate. This seems to be happenstance, since changes in 



the individual contributions are larger than the change 
in the total. The main changes occurred because we feel 
use of a larger magnetic polarizability is justified and be- 
cause using a dispersive treatment throughout does not 
allow keeping the elastic non-pole contributions. 
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